Ab-initio GMR and current-induced torques in Au/Cr multilayers 
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We report on an ab-imtio study of giant magnetoresistance (GMR) and current-induced-torques 
(CITs) in Cr/Au multilayers that is based on non-equilibrium Green's functions and spin density 
functional theory. We find substantial GMR due primarily to a spin-dependent resonance centered 
at the Cr/Au interface and predict that the CITs are strong enough to switch the antiferromagnetic 
order parameter at current-densities ~ 100 times smaller than typical ferromagnetic metal circuit 
switching densities. 
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I. INTRODUCTION 

Magnetic metals are often well described using the 
effective mean-field description provided by the Kohn- 
Sham equations of spin-density functional theory. In 
this description, the Kohn-Sham quasiparticles experi- 
ence exchange-correlation potentials with a spin depen- 
dence that is comparable in strength to band widths and 
other characteristic electronic energy scales. The spin- 
dependent part of the Kohn-Sham quasiparticle poten- 
tial acts like an effective magnetic field that is locally 
aligned with the electron spin-density. Because of these 
strong spin-dependent potentials, the resistance of a fer- 
romagnetic metal circuit will change substantially when 
the magnetization orientation in any part of the circuit 
is altered, an effect known as giant magnetoresistance 
(GMR). Conversely, transport currents can destabalize 
magnetization configurations that are metastable in the 
absence of a current and change collective magnetiza- 
tion dynamics. In the case of circuits containing fer- 
romagnetic elements the influence of transport currents 
on the magnetization can be understood as following 
from the conservation of total spin angular momentum; 
the torques that reorient quasiparticle spins as they tra- 
verse a non-collinear magnetic circuit are accompanied by 
current-induced reaction torques (CITs) that act on the 
magnetic condensate. This type of phenomenon is not 
by any means limited to magnetic systems. For exam- 
ple, there have been recent studies of the interaction be- 
tween transport and charge density waves (CDW) , which 
find that the CDW order parameter configuration influ- 
ences transport and conversely that transport can alter 
the CDW,^ 

GMR and CITs have been extensively studied, both 
experimentally and theoretically, in ferromagnetic metal 
circuits lii^ In the case of ferromagnets current-induced 
torques are generally referred to as spin-transfer torques, 
since they represent a transfer of conserved total spin an- 
gular momentum between the magnetic condensate and 



quasiparticle currents. The spin-transfer argument is 
very general but is perhaps a little vague in that it does 
not always specify what is meant by the magnetic con- 
densate, {the magnetization). (Some work does appeal 
explicity to a s-d picture of transition metal magnetism 
with d-electron local moments exchange coupled to itin- 
erant s-electron bands.) 

Several of us have recently proposed^ that CITs can 
also act on the order parameter of antiferromagnetic met- 
als. Our proposal follows from a microscopic picture^ of 
CITs that applies to a metal with any kind of magnetic 
order, including antiferromagnetic order, and suggests 
that CITs are a universal phenomenon in magnetic met- 
als. In this picture CITs arise from the dependence on 
bias voltage of the relationship between the steady-state 
Kohn-Sham quasiparticle density matrix and the Kohn- 
Sham single-particle Hamiltonian. This picture in effect 
explicitly identifies the energy region between the chemi- 
cal potentials of source and drain, the transport window, 
with the transport electrons and the energy region be- 
low this with the magnetic condensate. Given the Kohn- 
Sham Hamiltonian, the condensate contribution to the 
density matrix can be constructed^^ (for slow condensate 
dynamics) by solving the time-independent Schroedinger 
equation for electrons in a system of interest including if 
appropriate their coupling to electrodes in equilibrium, 
whereas the transport contribution to the density ma- 
trix is constructed by solving the Schroedinger equation 
with the scattering-theory boundary condition that the 
electrons be incident from the source. All this is ac- 
complished conveniently using a non-equilibrium Greens 
function theoryii^ The change in condensate dynamics, 
or the current induced torques, follow from the differ- 
ence between the Kohn-Sham Hamiltonian constructed 
from the equilibrium density-matrix and the Kohn-Sham 
Hamiltonian constructed from the density-matrix in the 
presence of a transport voltage. For magnetic conden- 
sate dynamicSfi^ what is relevant is changes in the spin- 
dependent parts of the Hamiltonian and in-particular 



changes in the direction of the exchange-correlation ef- 
fective magnetic fields on each site. 

Several of us have recently used this approach to ex- 
plore the possibility of CITs in circuits containing anti- 
ferromagnetic metals'^ by studying a simple one-band, 
Hubbard-interaction, toy model. For this model, we 
found CITs that drive the antiferromagnetic order pa- 
rameter which are, remarkably, proportional to film 
thickness, provided that inelastic scattering^ is ignored. 
(In ferromagnets spin-transfer torques saturate at a finite 
value for large ferromagnetic film thicknesses.) In this 
paper, we evaluate CITs for a potentially realistic mag- 
netoelectronic system using an ab-initio non-equilibrium 
Green's function formalism.— This calculation convinc- 
ingly demonstrates that GMR and CITs do occur in 
circuits containing only antiferromagnetic and paramag- 
netic elements, and that they can be large. The ab-initio 
approach accounts for all electronic structure details of 
the materials and the interfacial bonding between ma- 
terials, while the NEGF formalism enables the calcula- 
tion of finite bias properties, such as current, and non- 
equilibrium spin densities. The specific calculations we 
report on in this paper were performed on a system with 
antiferromagnetic (100) growth direction bcc Chromium 
(Cr) leads separated by a fee Gold (Au) spacer. The 
Au(100)/Cr(100) multilayer system we consider here ap- 
pears to be attractive as a model system for antifer- 
romagnetic metal spintronics. Au/Cr multilayers and 
(100) growth direction epitaxy were studied some time 
ago both experimentally^ and theoretically,^ motivated 
in part by superconductivity that can occur in disor- 
dered Au/Cr films. We find that both GMR and CITs 
do occur in this ferromagnet-free magnetoelectronic cir- 
cuit, as anticipated by previous work4 Our calculations 
also identify new physics not anticipated in the early toy- 
model study. The new features are associated with spin- 
polarized interface resonance at the Au/Cr interface and 
with the presence of more than one propagating Cr chan- 
nel at the Fermi energy for some transverse wavevectors. 
The paper is organized as follows: in Sec. |TT1 we briefly 
review our calculation method, while in Sec. Illli we give 
specific details of the system under consideration. We 
explain our results in Sec. |TV] and finally, in Sec. |V] we 
present our conclusions and suggest some possible direc- 
tions for future research directed toward exploring the 
potential of spintronics in Au/Cr and related systems. 



II. METHOD 

To calculate equilibrium and non-equilibrium prop- 
erties of the system, we employ the non-equilibrium 
Greens' functions (NEGF) formalism within the density 
functional theory framework. The central object in this 
formalism is the equal time lesser Green's function, de- 
fined as G'^j,{t,t) = i{cj,{t)cj{t)). Here the labels 
refer to the single-particle basis set in which the non- 
equilibrium density matrix is evaluated. In the NEGF 



formalism one of the labels must specify an atomic site, 
and the other labels specify included spin and the or- 
bitals included at each site. These calculations normally 
assume transverse periodicity, either on an atomic scale 
as in this calculation, or in the context of a transverse 
super-cell model for transport through a system with a 
finite cross section. The single-particle basis therefore 
also includes a transverse wavevector label. The NEGF 
formalism implementation employed here has been de- 
scribed in detail in previous work.^^ In the present work 
we assume that the time-dependence of the mean-field 
Hamiltonian, which is associated with slow condensate 
dynamics when present, can always be ignored in us- 
ing the formalism to calculate a steady state value of 
G^j,{t,t). This steady state Greens function specifies 
the non-equilibrium density-matrix p of the system. 

Once p is determined for a particular magnetiza- 
tion orientation configuration, the GMR is obtained 
by finding the spin-dependent transmission: = 
Tr[(Im(E2)G"^Im(E^)G°)^]. Here C^"' is the retarded 
(advanced) Green's function for the device, while S^j-^-j 
is the self energy which accounts for the semi-infinite left 
(right) lead. 

This is the first ab-initio evaluation of either GMR or 
CITs in a magnetoelectronic circuit that does not con- 
tain ferromagnetic elements. CITs in antiferromagnetic 
systems cannot be estimated on the basis of total spin 
conservation considerations because the antiferromagnet 
order does not carry a total spin. The CITs must be 
computed microscopically by finding the change in the 
local spin-dependent exchange-correlation potential on 
each site due to the presence of non-equilibrium, cur- 
rent carrying electrons, and summing over all sites assum- 
ing rigid antiferromagnetic order. Rigidity is maintained 
because the torques associated with short-distance-scale 
changes in relative spin orientation are very large com- 
pared to CITs. The CITs are significant experimentally 
because they compete only with much smaller anisotropy 
torques. This picture of current induced torques, and its 
implementation in the NEGF formalism, is spelled out in 
more detail in Ref. fl6|. 

The contribution of atom i to the CIT (per current) is: 

In the above a refers to an orbital, so that the sum is over 
all orbitals a of atom i. q, • r/2 is the spin-dependent 
part of the exchange-correlation potential for orbital a 
on atomic site i; its magnitude is in effect the size of 
the spin-splitting for this orbital and site. Here r is the 
vector of Pauli spin matrices. (The presence of more than 
one orbital on a site and the orbital dependence of the 
spin-splitting fields they produce is one of the important 
differences between more realistic models of transition 
metal magnetism and simple generic 1-band toy models 
of ferromagnetism.) The orbital and site-dependent non- 
equilibrium spin density which appears in Eq.([T]) is given 
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by 

=Tr[pi^/V]. J = (^,«) (2) 

with the trace in spin space. Since the system we consider 
is metaUic, we assume that the non-equihbrium quanti- 
ties in which we are interested can be evaluated in hnear 
response; specifically for the non-equilibrium density ma- 
trix: 

Ptr = G'■Im(EDG^ (3) 

with all quantities evaluated at the Fermi energy. The 
torque which acts on the staggered antiferromagnetic or- 
der parameter is the corresponding sum of the torques 
on individual atoms. Since the directions of the spin- 
dependent exchange-correlation fields alternate from site 
to site, current induced transverse spin-densities of the 
same sign give total torque contributions of the same 
sign>^ 

III. CR/AU/CR MULTILAYERS 

We study a circuit with semi-infinite antiferromagetic 
Cr leads and a Au spacer. Cr has a bcc lattice structure 
with lattice constant 2.88 A, while Au is fee with lattice 
constant of 4.08 A. The interface between these materi- 
als has a fortuitous lattice matching of two-dimensional 
square nets when they are grown epitaxially along the 
[001] direction and the Au lattice is rotated by 45° around 
the growth direction. In this configuration the bulk 
square net lattice constants differ by less than .2%. The 
lattice matching strains for Au on a [001] Cr substrate 
are therefore quite small. 

The antiferromagnetic state of bulk Cr has been stud- 
ied extensively ."iZ, The origin of antiferromagnetism in Cr 
is nesting between electron jack and hole neck pieces of 
the paramagnetic Fermi surf ace. J^J^ The nesting vector 
Q defines the spin density wave (SDW) period which 
is nearly commensurate with the lattice with Qaj^-n — 
.95, where a is the lattice constant. In thin film and 
multilayer structures, Cr can exhibit paramagnetism as 
well as commensurate and incommensurate SDW states, 
depending on the film thickness and on the adjacent 
materials There is evidence of antiferromag- 
netism in Cr thin films grown on Au substrates for cov- 
erages greater than 12 ML.^° Density functional theory, 
within LSDA, has been previously used to study bulk 
Cr in a commensurate SDW state^^ , and can predict the 
magnitude of the exchange splitting, the magnetic mo- 
ment, and the lattice constant. For this study, we restrict 
our attention to Cr with a commensurate spin density 
wave structure in which the body center spins and cube 
corner spins have opposite orientations. This magnetic 
structure is metastable in the absence of a bias voltage 
for our thin film structures in the local spin-density ap- 
proximation. All interfaces are perpendicular to the [001] 
direction. 
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FIG. 1: GMR as a function of spacer thickness. There is a 
sizable GMR for all spacer thicknesses. The inset shows the 
geometry for the four layer spacer configuration - up and down 
spins are colored red and blue, respectively. The configuration 
above illustrates an anti-parallel configuration. 

For the GMR calculations discussed below we have 
used a double-zeta with polarization basis set for both Cr 
and Au, and have found excellent agreement with bulk 
band structures and density of states calculations. For 
the calculation of current induced torques, we have used a 
single-zeta with polarization basis set, which still retains 
good accuracy for bulk properties. For the calculation 
of the self-consistent equilibrium density matrix, contri- 
butions from 900 fc-points within the 2-d Brillouin zone 
have been summed. For the calculation of the conduc- 
tance and current induced torques, 25600 /c-points have 
been summed. 

As a matter of convention, we define parallel (P) and 
anti-parallel (AP) alignment below in terms of the align- 
ment of the Cr spins in the two layers on opposite sides 
of the Au spacer. The inset of Fig. 1 shows the geometry 
and spin structure of the systems considered. 

IV. RESULTS 

Fig 1. illustrates the dependence of the GMR ratio 
(defined as (/p — Iap)IIap) on spacer thickness. In the 
limit of no spacer, the GMR ratio is negative, imply- 
ing larger conductances for an antiparallel configuration. 
This property is anticipated since the AP case at zero 
spacer thickness corresponds simply to ballistic conduc- 
tion through bulk antiferromagnetic Cr with unit trans- 
mission coefficient for all channels, whereas the parallel 
configuration implies a kink in the Cr antiferromagnetic 
order parameter configuration which reduces the trans- 
mission. For all nonzero spacers we have studied, we find 
that the GMR ratio is positive. The nonzero GMR for 
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FIG. 2: Layer and spin resolved density of states at the Fermi 
energy for Cr-Au single interface system. Layers 1-12 are 
Cr, and 13-24 are Au. The pronounced enhancement of the 
density-of-states in layer 12 reflects the interface resonance. 
The spin-polarization of the interface resonance is consistent 
with the fact that Tj > in the transport calculation. 



antiferromagnetic systems is perhaps surprising at first 
sight; for example a simple Julliere type two-channel con- 
ductor model, in which MR is due to spin-dependent con- 
ductance in the bulk, would predict that the GMR ratio 
is zero for antiferromagnetic systems. For antiferromag- 
nets GMR is, in this sense, purely an interface effect; for 
ferromagnets GMR is only partly (but often mainly) an 
interface effect. 

In the case of the previous toy model antiferromagnetic 
calculations^., GMR was due to phase coherent multiple 
scattering between two antiferromagnets. These effects 
are partially mitigated'' at elevated temperatures by in- 
elastic scattering which breaks phase coherence. The 
present ab-initio calculations reveal a new contribution 
to antiferromagnetic GMR, explained below, which does 
not rely on phase coherence. The property that realistic 
antiferromagnets have GMR effects that are not depen- 
dent on phase coherence is encouraging from the point 
of view of potential applications, since it suggests larger 
robustness at elevated temperatures. 

In order to identify the dominant GMR mechanism of 
Au/Cr multilayer systems, we have performed a separate 
NEGF calculation for a single interface between semi- 
infinite bulk Cr, and semi-infinite bulk Au. The result is 
that there is a spin dependent conductance with magni- 
tude (/| - Il)/{I^ + Ii) = -2.10%. The current is spin- 
polarized in the direction opposite to the top layer of the 
antiferromagnet. For Cr/Au/Cr multilayers, this spin- 
filtering implies that the conductance is maximum when 
the facing layers of the antiferromagnet have the same 
spin-orientation, i.e. the P configuration has a higher 
conductance, even without any non-local coherence ef- 



FIG. 3: Number of propagating states in the [001] direction 
at the Fermi energy for bulk antiferromagnetic Cr. 

fects. This effect is absent'^ in the single-band models 
that we studied earlier. 

To explore the origin of the spin-dependent interface 
resistance in greater detail, we have examined the layer 
and resolved local density of states of the single- 
interface calculation (here fcy refers to the transverse mo- 
mentum label). Fig. 2 shows the layer-resolved results. 
(In the calculation, the first 12 layers on each side of 
the interface are allowed to differ from the bulk. In Fig. 
2 layers 1-12 are the Cr layers in the scattering region 
near the interface while layers 13-24 are the Au layers 
in the scattering region.) We see that there is a pro- 
nounced interface resonace on the last Cr layer; this is a 
consequence of the difference between the Fermi surface 
topologies of Cr and Au. Moreover, this state is spin po- 
larized, with direction opposite to that of the bulk local 
moment. Fig. 3 shows the number of propagating chan- 
nels in the Brillouin zone for Cr, which demonstrates that 
the Fermi surface of Cr differs strongly from the nearly 
spherical Fermi surface of Au. In particular, there are 
large regions in the Brillouin zone of Cr in which there is 
no propagating state, whereas Au has propogating states 
across all of the central region of the transverse Brillouin 
zone. Fig. 4 shows the transverse wavevector-resolved 
Fermi-level local density of states for layers 8, 10, 12 (the 
last Cr layer), and 16. Layer 8 is typical of bulk Cr, while 
layer 16 is typical of bulk Au. Layer 12, however, shows 
features of both materials; in particular, populations of 
states within the region of the Cr Brillouin zone with no 
propagating modes are responsible for the localized in- 
terface resonance. Fig. 5 shows the total local-density of 
states as a function of energy for layers near the inter- 
face. The rapid relaxation toward bulk values away from 
the the interface is apparent. The interface layer has a 
highly distorted density of states function, a high density 
of states at the Fermi level, and a higher moment den- 
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sity which is responsible for a net ferromagnetic momen-bs 
contribution from the interface region. Apparently in- 
terruption of antiferromagnetic order both narrows the 
majority-spin bands and lowers the energy of minority 
spins in this interface layer. Hopping of down-spins from 
the sub-interface layer on which they are the majority to 
the spacer layer is enhanced by the minority spin inter- 
face resonance. 

The enhanced moment density in the interface layer 
is accompanied by more attractive spin-up potentials on 
this layer and spin-dependent bonding across the Cr/Au 
interface. The effective hopping matrix elements across 
the interface have a spin dependent contribution that 
is about 1% of their total values. In order to deter- 
mine whether it is spin-dependent hopping or resonances 
related to spin-dependent site energies we have sym- 
metrized hopping of the interface to remove its spin de- 
pendence, and re-calculated the conductance. We find 
the same value for the polarization, indicating that it is 
the interface resonance that is largely responsible for the 
polarization. The fact that the GMR is due to inter- 
face resonances, rather than to phase-coherent multiple- 
scattering across the spacer layer, suggests that the effect 
will be robust at elevated temperatures. 

As mentioned earlier, the antiferromagnet/normal in- 
terface resistance is not spin-dependent in the toy model 
systems previously studied^' The key property of the toy 
model which leads to this spin-independent interface re- 
sistance is that each antiferromagnetic unit cell is invari- 
ant under a combination of space and spin inversion. The 
ab-initio mean-field Hamiltonian does not have this prop- 
erty. A spin-dependent resistance will occur in the toy 
model when either the hopping amplitude from spacer to 
the top antiferromagnetic layer is made spin-dependent 
or the exchange-splitting in the top layer is shifted from 
its bulk value. 

We have evaluated the current induced torques for a 
system with a four Au monolayer spacer. The angle be- 
tween the staggered moments of the Cr leads was initial- 
ized to 90°: the staggered magnetization is along the z 
direction in the Cr layers to the left of the spacer and 
along the x direction in the Cr layers to the right. A 
self-consistent non-collinear solution to the Kohn-Sham 
equations was obtained with this configuration. The re- 
sulting layer resolved torques, evaluated as described in 
Sec. ini are plotted in Fig. 6. We find strong torques 
peaked in the first few Cr layers, in contrast to the toy 
model case in which the torques were constant in mag- 
nitude and alternated in direction from layer to layer. 

To shed light on the origin of the new antiferromag- 
net current-induccd-torque physics revealed by these ab- 
initio calculations we focus on differences between the 
toy-model case, in which analytic calculations are pos- 
sible, and the realistic Au/Cr case. Because the equi- 
librium torques that impose the antiferromagnetic order 
will always be much stronger than the current-induced 
torques, we are mainly interested in the sum of the layer 
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FIG. 4: Transverse-momentum resolved Fermi-energy 
density-of-states 8, 10, 12, and 16 of the single interface Cr- 
Au system. Layer 12 is the Cr layer closest to the interface. 
Layer 8 shows bulk Cr characteristics (compare to Fig. 3), 
while layer 16 shows bulk Au characteristics. 



resolved torques which drives the antiferromagnetic order 
parameter and therefore competes only with anisotropy 
torques. The perfectly staggered torque obtained in 
the toy model case arises from an out-of-plane current- 
induced spin density that is spatially constant within 
each antiferromagnet. (Here in-plane refers to the plane 
spanned by the orientations of the two antiferromagnetic 
layers and out-of-plane refers to the perpendicular di- 
rection, the y-direction in our case.) The constant out- 
of-plane current-induced spin density in the toy model 
can be partially explained by the fact that Bloch wave 
vectors of up and down spin states are not spin-split in 
antiferromagnets. It follows that a linear combination of 
transmitted up and down spins have a transverse spin 
density that is position-independent. (In contrast to the 
ferromagnetic case, in which the transverse spin density 
of a particular channel shows spatial oscillations with a 
period given by (fc| — /c|)^^).^^ To see where this physics 
breaks down in our calculation, consider Fig. 3, which 
shows that a particular transverse channel has 1-4 possi- 
ble values of ■ For those channels with a single value 
at the Fermi energy, we find that the contribution to the 
transverse spin density is spatially constant. Evidently 
the toy model does a good job of describing this type of 
transverse channel, suggesting that our earlier conclusion 
that there is a bulk contribution to the staggered spin- 
torque in an antiferromagnet does have general validity. 
The present calculations emphasize, however, that there 
is also an interface contribution coming dominantly from 
channels with more than one kz value at the Fermi en- 
ergy. In this case the transmitted wave function is a lin- 
ear combination of states with different fc^ values. These 
states interfere with each other to produce an oscillat- 



5 



ing transverse spin density. Summing over many chan- 
nels with different osciUation periods of the transverse 
spin density leads to a rapid decay of the transverse spin 
density, exactly as in the case for ferromagnets.^- A ma- 
terial with a simpler, single valued Fermi surface would 
not have this interface contribution to the total staggered 
torque. The complex Fermi surface is necessary to stabi- 
lize the AFM order via nesting in the first place, however. 
In our case, the interface torque dominates the total stag- 
gered torque. It is nevertheless substantial, totalling .049 

To estimate the critical current for switching the anti- 
ferromagnetic order parameter, we take the anisotropy 
of bulk Cr with spins pointing along the n direction 
as: E{n) = ifi(z • n)^ -|- 4is:2(x • n)2(y • n)2, wherefi^ 
Ki = 10'^ J m^'^, and ^2 = 10 J m"-^, and take the 
magnetic damping parameter to be a — 0.1.— Here Ki 
is positive for T > Tsf = 123. 5if, and Q, the spin 
density wave vector, is taken to be in the z-direction. 
Near the fixed point, n = x, the damping torque per 
area is then F = aj{Ki + AK2)t, where t is the thick- 
ness of the layer. (Note that antiferromagnets possess 
no demagnetizing field, so that the anisotropy does not 
include shape anisotropoy and is due only to magne- 
tocrystalline anisotropy.) The current required to have 
the current-induced-torque overcome damping is there- 
fore M9{^t/e)r « 6.3-lO^H{A/m^). (In this we assume 
that the transfer torque efficiency calculated at 9 = 90° 
is the same as that for small or large angle.) Typical 
values for critical current densities in ferromagnets are 
up to 100 times larger, primarily because of the large de- 
magnetizing fields present in ferromagnets. We therefore 
expect that it will be possible to achieve current induced 
switching in a circuit containing only antiferromagnetic 
elements. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have performed a first-principles cal- 
culation of the non-equilibrium properties of Au/Cr mul- 
tilayers. We have found that a system composed of Cr 
leads with an Au spacer possesses a robust GMR effect, 
which is due primarily to a spin polarized interface state 
at the Au/Cr interface, and a strong current induced 
torque which corresponds to a switching critical current 
more than an order of magnitude smaller than that of 
ferromagnets. These ab-initio calculations demonstrate 
that current-induced-torques in antiferromagnets have a 
bulk contribution from channels with only one band at 
the Fermi energy and an interface contribution, which 
is dominant in our case, from channels with more than 
one band at the Fermi energy. Our calculations demon- 



strate that all the basic effects of metal spintronics occur 
in circuits with only paramagnetic and antiferromagnetic 
elements. 

The robust antiferromagnetic spintronics effects stud- 
ied here occur in circuits in which current ffows per- 
pendicular to planes containing perfectly uncompensated 
spins. We believe that substantial effects will also occur 
in realistic experimental systems in which all layers are 
partially compensated. Because the order in an antiferro- 
magnet is staggered, structural disorder (interface rough- 
ness for example) will induce magnetic order, "^^ Although 
we have not studied these effects in any depth to date, it 
seems likely that observable GMR and current-induced 
order-parameter and resistance changes will survive as 
long as the length scales associated with spatial disor- 
der are larger than the spacer layer thickness. For Cr 
on Fe whiskers the lateral correlation length is the ter- 
race width, which can be easily larger than typical spacer 
widths if the interface is fiat3^i2ii2^ 

It may be that initial studies of spin-torques in antifer- 
romagnetic systems will be easier to analyze in systems 
in which the antiferromagnets are exchange coupled to 
ferromangets, whose orientation can be manipulated by 
external magnetic fields. There is, for example, already 
evidence^ that current-induced torques act on the an- 
tiferromagnetic layer in spin-valve structures. Systems 
containing Cr antiferromagnetic films that are exchange 
coupled to Fe whiskers might also be attractive to study 
these kinds of current-induced torques since it is already 
known that weakly compensated Cr layers can be ob- 
tained relatively easily.— i^iii^ It may also be advanta- 
geous to consider Mn doped Cr, as it is known to have 
a higher Neel temperature than pure Cr, and forms a 
commensurate spin density wavci^ Although the ma- 
terials challenges presented by antiferromagnetic metal 
spintronics are even stronger than those presented by fer- 
romagnetic spintronics, we believe that the subject will 
prove interesting from both basic physics and potential 
application points of view. 
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FIG. 5: Density-of-states for layers 10-14 of the single interface Cr-Au system. The DOS relaxes to its bulk shape a couple of 
layers away from the interface. Layer 12 is the interface Cr layer, end the Fermi energy is 0. 
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